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It is demonstrated that low frequency shear modes in a strongly coupled, inhomogeneous, dusty 
plasma can grow on account of an instability involving the dynamical charge fluctuations of the 
dust grains. The instability is driven by the gradient of the equilibrium dust charge density and is 
associated with the finite charging time of the dust grains. The present calculations, carried out in 
the generalized hydrodynamic viscoelastic formalism, also bring out important modifications in the 
■ threshold and growth rate of the instability due to collective effects associated with coupling to the 
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compressional mode. 
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I. INTRODUCTION 



Dusty plasmas are of great interest because of their possible applications to a number of fields of contemporary 
research such as plasma astrophysics of interplanetary and interstellar matter, fusion research, plasmas used for 
semiconductor etching, arc plasmas used to manufacture fine metal and ceramic powders, plasmas in simple flames 
, etc. Q. It is now widely recognized that the dust component in these plasmas is often in the strongly coupled coulomb 
regime with the parameter, T ~ (Z^e) 2 /Tjd, typically taking values much greater than unity ((— Z^e) is the charge 

— 1/3 

on the dust particle, d ~ n d is the interparticle distance and Td is the temperature of the dust component). This 
leads to many novel physical effects such as the formation of dust plasma crystals modified dispersion of the 
compressional waves [pLW, the existence of the transverse shear waves M etc. Many of these novel features have now 
been verified by experiments and computer simulations Q . 

7-H \ Recently, an experiment on the self-excitation of the vertical motion of the dust particles trapped in a plasma 
J> ■ sheath boundary, has been reported |6| . The physics of this excitation is related to charging of the dust particles by 
the inflow of ambient plasma currents in the inhomogeneous plasma sheath and the delay resulting because of the 
finite time required by the charging process to bring the dust charge to its ambient steady state value. In this paper, 
we demonstrate that the same physical mechanism can be used for the excitation of the transverse shear modes in 
an inhomogeneous strongly coupled dusty plasma. Using a generalized hydrodynamic viscoelastic formalism jTj] to 
' describe the strongly coupled dusty plasma and incorporating the novel feature of time variation of the dust charge 
through a charge dynamics equation we have derived a general dispersion relation for low frequency shear and 
compressional modes in the plasma. We find that in a plasma with finite gradients of the equilibrium dust charge 
density, the two modes are coupled and we show that the shear mode is driven unstable if certain threshold values 

" ^ ■ are exceeded. 

& Our paper is organized as follows. In the next section we briefly discuss the equilibrium of an inhomogeneous dusty 

plasma that is confined against gravity by the electric field of a plasma sheath. In such a configuration dust particles 
of varying sizes and charges arrange themselves in horizontal layers at different heights to form a nonuniform cloud 
P,flO[ . In section 3 we carry out a linear stability analysis of such an equilibrium in the framework of the generalized 
hydrodynamic equations. The dispersion relation of the coupled shear wave and compression wave is solved analyt- 
ically (in simple limits) as well as numerically in section 4. The physical mechanism of the shear wave instability is 
also discussed and the modifications in the threshold and growth rate brought about by the coupling to compressive 
waves are elucidated. Section 5 is devoted to a summary and discussion of the principal results. 
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II. DUST CLOUD EQUILIBRIUM 



We consider an inhomogeneous sheath equilibrium in which the dust particles are suspended with electric field 
forces balancing the gravitational force on the particle and in which the dust charge (— Z^e) and dust size r& are both 
functions of the vertical distance z. Then the force balance equation gives, 

Z d (z)eE (z) = ^irr d (z) 3 pg, (1) 

where, p , g , Eq refer respectively to the dust mass density, gravitational acceleration and the sheath electric 
field. For particle sizes of the order of a few microns, other forces acting on the particle (such as the drag and viscous 
forces) are about an order of magnitude smaller than the electric and gravitational forces and can be neglected for 
the equilibrium calculation fT(f | . Note that for dust particles of a uniform size (monodispersive size distribution) the 
above equilibrium can only be attained at one vertical point leading to a monolayer of dust. A dispersion in sizes leads 
to a large number of layers resulting in a nonuniform dust cloud with a gradient in the equilibrium charge (—Z^e) 
and the dust size r d . The electric field E is determined by the sheath condition, 

dE$ 

— — = -47re (n e - n t + ZdUd) (2) 

where n e ^ y d are the local electron, ion and dust densities respectively. The charge (— i^e) on a dust particle in 
the sheath region is given by (— Z#e) — Cd{4>f ~ 4 1 ) where Cd is the capacitance, <f>f is the floating potential at the 
surface of the dust particle and <f> is the bulk plasma potential. For a spherical dust particle Cd = rd, and the floating 
potential can be determined by the steady state condition from the dust charging equation, namely, || 

h + h = (3) 
where the electron and ion currents impinging on the dust particle are given by |l[ 



I e = -nr^e ( I n e exp 



8kT e 



1/2 



(4a) 
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I, = 7T7\,e 



1/2 

(4b) 
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Here T e and Ti are the electron and ion temperatures, mi is the ion mass and Vid is the mean drifting velocity of the 
ions in the electric field of the sheath (it is assumed to be the ion sound velocity at the sheath edge). We also assume 
that the dust particles have much smaller thermal velocities than the electrons and ions. 

Equations (|l| - ||) selfconsistently determine the equilibrium of the dust cloud. Such clouds have been experimentally 
observed in a number of experiments p|Jl0[|. In |lCf| theoretical modeling along the lines discussed above, agree very 
well with the experimental observations of clouds formed with polydispersive particle size distribution of dust particles 
trapped in the plasma sheath region. A typical equilibrium variation of the dust particle size with the vertical distance, 
when the Child Langmuir law holds for the plasma sheath potential, is given as JlC[ ], 

rd = ( Wf-*)\ V2 ( 8«en a C B \ 1/3 {§ _ 1/8 (g) 



^pg J V M 

where n s , C s are the plasma density and the ion sound velocity, S is the thickness of the sheath and 
Pi = (eXi-n/rrii) 1 / 2 with \i- n representing the mean free path of ions colliding with the background neutrals. Using 
(|5|) we can obtain the corresponding z variation for Zd- 

As discussed in detail in (h]] , this dust cloud equilibrium is confined to the plasma sheath boundary region in the 
potential well created from the upward electrostatic and downward gravitational forces. Note that the force balance 
equation ([!]) does not prevent the particles from oscillating about their mean positions especially if they have signifi- 
cant kinetic energy or temperature. However their mean positions are at various vertical distances and the mean Zd 
is a function of z. This is reminiscent of particle gyrations in a magnetic field. If we consider wave motions in which 
dust oscillation excursions are much smaller than wavelengths, we can use a fluid theory to analyze such behaviour. 
In the next section, we adopt this view point and carry out a linear stability analysis of the equilibrium discussed 
above to low frequency wave perturbations. 
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III. LINEAR STABILITY ANALYSIS 



For low frequency perturbations in the regime kvthd << <*> << kv t he,kvthi, where v t hd, Vthe and Vthi are the 
thermal velocities of the dust, electron and ion components respectively, the electron and ion responses obey the 
Boltzmann law which can be simply obtained from an ordinary hydrodynamic representation. The dust component 
on the other hand can be in the strongly coupled regime for which a proper description is provided by the generalized 
viscoelastic formalism. Using such a description a general dispersion relation for low frequency waves (with typical 
wavelengths longer than any lattice spacings) was obtained in B for longitudinal sound waves and transverse shear 
waves. The shear modes exist in a strongly coupled dusty plasma because of elasticity effects introduced by strong 
correlations Our objective in this work is to look for the effect of dust charge dynamics on these shear modes in 
the strongly coupled regime. As demonstrated in our earlier work the coupling of the low frequency shear modes 
to transverse electromagnetic perturbations is finite but negligibly small; we ignore this coupling here. However, 
introduction of the dust charge dynamics in the inhomogeneous plasma leads to a coupling of the low frequency shear 
and compressional modes; thus the space charge dynamics and quasineutrality condition play an important role in 
describing the perturbations. The basic equations for the dust fluid we work with, are the continuity equation, 



^-5n d + n d0 V ■ Su d + ^-(5u d • V)M = 0, 
at M 



(6) 



the equation of motion, 



X/SP , Z d e 



Mn d0 M 



SE 



5Zd 
M 



■eE 



Mn d0 



riy 2 6u d + (C + o)v(V ' <M 



and the equation of state, (dP/dn)? — MC%, given in terms of the compressibility, /j, d , as 



= ±(—) = 1 | " (F) | Tdu{r) 



T d \dnJT * ' 3 ' 9 dT 
with the excess internal energy of the system given by the fitting formula fll|| 

u(T) = -0.89r + 0.95r 1/4 + 0.19r~ 1/4 - 0.81. 



(7) 



(8) 



(9) 



In the above, M is the dust mass, v is the dust-neutral collision frequency, 5u d , Sn d and SZ d are the perturbations 
in the dust velocity, number density and dust charge, 5P, SE are the pressure and electric field perturbations, n d Q 
and Z d are the equilibrium number density and charge for the dust and Eq is the unperturbed electric field, rj and £ 
refer to the coefficients of the shear and bulk viscosities and r m is the viscoelastic relaxation time. Note that in the 
continuity equation we have a contribution from the equilibrium inhomogeneity in the dust mass distribution (arising 
from the size dispersion of the particles). This term as we shall see later modifies the real frequency of the shear 
waves. 

These equations are supplemented with the dynamical equation for the dust charge perturbations which, for per- 
turbations with phase velocity much smaller than the electron and ion thermal velocities, is given as |q| 



— (5Z d ) + Su d ■ \/Z d + i] c bZ d = —- 

at e \ n i0 n e0 



(10) 



where, rj c = ^e|/ e o|/C^ (l/T e + l/wo \ is the inverse of charging time of dust grains and wq — Ti — e(cf>f — </>)o. Note 

that the second term on the left hand side of eq.(|l(]) arises because of the inhomogeneity of the mean charge on 
the dust particles; as shall be shown later, this is the critical term responsible for the instability. It is also obvious 
that the dust charge variation in space will lead to shielding by electrons and ions with the associated coupling 
of the perturbation to dust compressional modes. We must thus extend the above set of equations to include the 
quasi-neutrality condition, 



5n e + Z d Sn d + n d0 SZ d - 5m ~ 0, 



(11) 



3 



and the equation describing the electron and ion density perturbations in terms of the potential, as 



6th 
n i0 



(12) 



These are the Boltzmann relations which arise whenever the perturbations satisfy u> « kvthe, kvthi- 

We shall next derive the dispersion relation for the low frequency mode. We may note that the typical time scale 
for the decay of the charge fluctuations for the dust can be very small [||, with r\ c >> u> and we shall work in that 
limit. We use the local approximation (wave lengths smaller than characteristic equilibrium scale lengths) and choose 
the propagation vector for the wave perturbation as k = [k, 0, 0), the perturbed dust velocity, Sud = (Su\, 0, Su 3 ) and 
the perturbation in the electric field as SE — —ik8(p(l, 0, 0). Using the continuity equation (||) and the equations ( |l(i| ) 
- (|12|), and after some simple algebra, one obtains the fluctuation in the dust charge and the potential as 



5Zd ^D\u) 5ul + \D + J^)D 



Sun 



(13a) 



eD 



— OUl H Zj ; — — 

ujJ eD\ d M(iu>) 



5u 3 , 



(13b) 



where, 



ai = 



a 3 - 
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1 1 \ M' 

n + T, Ji ndaZd 



D = Vc 



n e o , riio 



T, 



ndo 



7 i / n e o n-io 



T, 




T, 



(14) 



and the primes denote derivatives with respect to z the vertical direction. We then write down the longitudinal and 
transverse components of the dust momentum equation (i.e. of equation (^)), as 



(1 - iujT m ) (~iuj + v)8ux + ifc-r?~ — ~!nr(^0) 

Mn d o M 



1 



Mn d o 



rjik Sui 



(15a) 



8 Z 

(1 - iuT m ) (-iw + v)5u 3 + -j^j-eEo 



Mn d0 



■qk 2 8u 3 , 



(15b) 



where, rji = |r/ + C. In the limit uiT m >> 1, using equations (|l3|)- (^5[), we obtain the dispersion relation for the 
coupled shear-compressional mode, as 



2 ■ e ^o 

U + IUJV + IU1 02 

2 eJ?p aiZdZ' d ndo 
,j)k MD MB + 



eE 
MD 



d3 



^ DA 



k 



2 eE aiM' 
MD M 



K + cl) = 0, 



(16) 



where C 2 sh = (rj / Mn d0 T m ) , C\ a = {Z 2 d n d ^c/MD) and C 2 DA 



(rji/MndoTm)- In the above equation the 



expression in the first set of brackets represents the dispersion relation for the transverse shear wave, the second set 
of brackets contains the compressive mode dispersion relation and the final two terms denote the coupling between 
the two branches. We will now study the behaviour of the shear mode in the presence of the charge inhomogeneity 
and the coupling to the compressive mode. 
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IV. SHEAR WAVE INSTABILITY 



In the limit when the coupling to the compressive wave is weak, so that the last two terms in the dispersion relation 
( pf ) can be neglected, we can obtain the roots for the shear branch as, 



if , e£ \ , f. 2^2 eE If eE \ 2 1 1/2 
■2\ V+ MD a2 ) ± [ k C °»~MD a3 --A" + MD a2 ) \ ' 



In the absence of the inhomogeneities and the collision term, this is the basic shear wave described in 0. The 
collisional term introduces wave damping. The inhomogeneous terms introduce two important modifications. The 
term proportional to the mass (size) inhomogeneity contributes to the real part of the frequency whereas the charge 
inhomogeneity term can drive the wave unstable if E a 2 < (i.e., E Q' < 0) and the threshold condition v < IfH^I 
is satisfied. Physically, this instability arises because of delayed charging effect, the same physical mechanism which 
was used by Nunomura et al ra] to explain the observed instability of single particle vertical displacement in their 
sheath experiments. Specifically, the charge on the vertically oscillating dust particle in the shear wave propagating 
in the inhomogeneous plasma, is always different from the equilibrium value Zd because of the finite charging time 
rj^ 1 . This perturbation is of order SZd — Z' d 5u^/r] c and leads to an energy exchange between the shear wave and 
the ambient electric field at a rate SZdEoSu^ « \EpZ' d \\8u^\ 2 /ri c . When this energy gain by the shear wave exceeds 
the loss rate due to collisions rs ^y~\8u^\ 2 , we have an instability. This gives us the approximate threshold condition 
described above. If we express the dust neutral collision frequency, v in terms of the ambient neutral pressure as 

v = p(^p !L ) 1 ^ 2 ' ^ffi °ur threshold condition is functionally identical to that derived by Nunomura et al |J on the basis 
of physical ar gum ents. The only substantial difference is their use of exponential charging time which follows from 
our equation ( |l0| ) viz. SZd ~ {8u^Z' d / rjc)\\ ~ G x p(~ilct)]j since we have assumed the frequency of the shear mode 
lu << r] c , we use the asymptotic condition described above. 

We now demonstrate that for the collective shear mode being described here, the coupling to the compressional 
dust acoustic wave due to the last two terms in equation ( |l6| ) is very crucial; thus the above single particle results are 
strongly modified by the hydrodynamic treatment. A simple analytic result clearly demonstrating the modification is 
obtained by neglecting to 2 + iujv compared to k 2 Cj-, A in the second bracket of equation ([H]); this is reasonable when 
the wave-vector k is not too small. In this limit, the shear modes are described by the root 

eE / ZdZ' d ndo 



U> = — -I V + —— 1 "2 " "j 



2 V 1 MD V" ' " x MD C 2 DA 



1/2 



sh MD\ M C 2 DA ) 4V MD\ MD C 2 DA > 7 J v ; 

We thus note that the threshold condition and the growth rates are significantly modified by the inclusion of coupling 
to compressional waves. In order to quantitatively illustrate the effect of coupling terms, we now present a detailed 
numerical investigation of the dispersion relation equation (|l6|). It is generally the case that the bulk viscosity 
coefficient £ is negligible compared to the shear viscosity coefficient, 77, particularly in the one component plasma 
OCP) limit and so we shall drop it in our calculations. Further, the viscoelastic relaxation time, r m , is given as 



4// 



3n d0 7d(l - IdfJ-d + j5«) 



4TT (19) 



with 7^ as the adiabatic index and the compressibility, fid defined through (||). We assume the gradient of the 
cquilibriated dust charge to be of the form, Z' d = Z^/Lz, the mass gradient to be of the form M' = M/Lm where 
Lz ~ Lm = L is a few Debye lengths. In our computations, we choose L rj 5 times the Debye length, which is the 
typical order of magnitude as observed experimentally Jl0[ . For further computations, we introduce the dimensionless 
quantities, 

u = uj/uj p d; v^v/ujpd; k = kd; r m = T m u pd ; 

V = 77 — " &a = Cl/{ujldd 2 )\a = sh,d,da,DA, 

MridoLUpdd^ 

_ eE a 2 _ eE a 3 

_ aiZ d Z' d n m 1 _ eE M' - 2 - 2 1 

601 - a 2 MD co 2 pd d 2 > eil ~MD ai J¥ {Cd+Cda) ^f d > (20) 
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where Lu p d and d are the dust plasma frequency and the inter-grain distance respectively. The dispersion relation for 
the shear mode ( |l6| ) can then be written as 

uj 2 + iCjv + iujeo + ei — C 2 h k 2 

— iujk 2 eoeoi + ke n = 0, (21) 

Equation ( pl| ) has been solved numerically for the shear mode roots and some typical results are presented in figures 
(1) and (2). Figures (la) and (lb) display a comparison of the dispersion curve for the shear mode (u>r vs k and 7 vs k 
for fixed values of eo = —0.0008 and e\ = —0.05), with and without the inclusion of the coupling to the compressional 
mode. The various fixed parameter values corresponding to these curves are C 2 h = 0.02, C\, A = 0.4, v = 0.0004 and 
eoi = 0.3, en = —0.01 when the coupling is on. The choice of these numerical values for the dimcnsionlcss parameters 
P, k, eo, eoi, C s h and C'da has been guided by the magnitude of these quantities observed in some of the laboratory 
plasmas [^|,|l0|. It is seen from these plots that there is a substantial influence of the compressional mode coupling, 
described through the parameter, eoi, en, on the growth rate and the real frequency of the shear wave emphasizing 
the importance of the collective physics of coupling to the compressional mode. We next plot in figure (2) the gas 
pressure, p, versus n e Q profiles for various values of 7, the imaginary part of tu. Plotting the 7 = curve, we get a 
threshold relation between p and n e o, where we fix the other parameters as follows - dust radius, ra=2.5 microns, the 
inter-grain distance, g?=430 microns, T e = Tj ~ leV, kd = 1, and dust mass density, pd—2.5 gms/cm 3 . We see that 
the qualitative trend of the curve is similar to that observed in the single particle instability studies of [[| illustrating 
the commonality of the underlying physical mechanism. However it should be emphasized that the experiment in 
f| did not observe any collective excitations and their equilibrium consisted of a monolayer of equal sized particles. 
The equilibria of |],[^| are more appropriate for observing collective excitations of shear waves and our theoretical 
results can be usefully employed in such a situation. In Fig. (2) we have once again highlighted the significance of 
the coupling to the compressive wave, in this case for its effect on the threshold values, by displaying the uncoupled 
threshold and growth rate curves (dashed curves). Note that the influence of the coupling is to raise the threshold 
value at low values of n e o (i.e. a higher value of p is needed to excite the instability) whereas it reduces the threshold 
at the higher end of the n e o scale. The rest of the curves displayed in the figure (2) correspond to the various positive 
values of 7, which correspond to the situation where the shear mode is excited and saturates at some values. These 
figures are again qualitatively similar to the curves obtained in [|| for various saturation amplitudes. However a direct 
comparison is again not appropriate for the reason discussed above and also because our calculations are linear and 
cannot provide any quantitative results about nonlinearly saturated amplitudes. 



lu 2 + iuji> — C nA k 



V. CONCLUSION AND DISCUSSION 



To summarize, in this paper we have investigated the stability of a low frequency shear mode in an inhomogeneous 
dusty plasma in the strongly correlated regime. The equilibrium dust cloud has both an inhomogeneity in the dust 
charge distribution and in the dust mass distribution (arising from a distribution in the sizes of the dust particles). 
The shear mode in such a plasma undergoes two significant modifications. Its real frequency is shifted by a contri- 
bution from the mass inhomogeneity and the dust charge inhomogeneity can drive it unstable through the dynamics 
of dust charge fluctuations in a manner very similar to the instability of the vertical motion of single particles in 
a plasma sheath as observed in the recent experiment of Nunomura et al The finite charging time, r\~ 1 of the 
dust particles plays a critical role in the instability. We also show how collective effects due to coupling with the 
compressional modes strongly modify the threshold conditions for the instability as well as its growth rate and real 
frequency. Our calculations have been carried out in the hydrodynamic formalism including viscoelastic effects and we 
have neglected any kinetic effects. Our results are therefore strictly valid in the low frequency limit. Finite corrections 
arising from kinetic effects can occur at higher frequencies and wave numbers. This has recently been demonstrated 
for the compressive dust acoustic mode in a dusty plasma from a kinetic calculation based on the dynamic local field 
correction (DLFC) method fl2[| . Such corrections, if any, for the transverse shear mode has not yet been done and 
needs to be examined. 

Finally we would like to remark that the transverse dust shear mode which is a collective mode of the strongly 
coupled plasma regime has only been observed in computer simulations till now; its detailed experimental investigation 
is therefore of great current interest. Such waves can be excited in inhomogeneous dust clouds that have been obtained 
in the experiments carried out with varying grain sizes ||9|,|l0| . It would be of interest therefore to look for the wave 
features discussed in our model calculations in controlled propagation experiments on such equilibria. It is also 
apparent that free energy sources, such as ion beams, which may readily couple with the compressional waves may also 
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be useful for exciting the more interesting shear waves in the strongly coupled inhomogcncous plasma. Investigation 
of these and related effects are in progress. 
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FIG. 1. (a) The normalized real frequency and (b) the normalized imaginary frequency, vs. the normalized wave number for 
the shear mode with e = —0.0008, e i = 0.3, e\ = —0.05, en = —0.01 (solid curves). The dashed curves are for e i = en = 
and correspond to the uncoupled shear mode. 
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FIG. 2. The electron number density n e o (in units of 10 /cm ) is plotted as a function of the gas pressure, p (in mtorr) 
for various values of 7. For comparison, the accompanying dashed curves display the situation when the coupling to the 
compressional mode is neglected. 
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